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1 Introduction 

The perturbations of a black hole are governed by quasi-normal modes (QNMs). The latter are 
typically obtained by solving a wave equation for small fluctuations in the black hole background 
subject to the conditions that the flux be ingoing at the horizon and outgoing at asymptotic infinity. 
These boundary conditions in general lead to a discrete spectrum of complex frequencies whose 
imaginary part determines the decay time of the small fluctuations 

3w = - (1) 

T 
There is a vast literature on quasi-normal modes and I make no attempt to review it here. Instead, 
I concentrate on obtaining analytic expressions for QNMs of black hole perturbations in asymp- 
totically AdS space. One can rarely obtain analytic expressions in closed form. Instead, I discuss 
techniques which allow one to calculate the spectrum perturbatively in an asymptotic regime (high 
or low overtones). For high overtones, the frequencies at leading order are proportional to the ra- 
dius of the horizon. For low overtones, one in general obtains an additional frequency which is 
inversely proportional to the horizon radius. Thus for large black holes there is a gap between the 
lowest frequency and the rest of the spectrum of QNMs. I pay special attention to the lowest fre- 
quencies because they govern the behavior of the gauge theory fluid on the boundary according to 
the AdS/CFT correspondence. The latter may have experimental consequencies pertaining to the 
formation of the quark-gluon plasma in heavy ion collisions. Moreover, I discuss phase transitions 
to hairy black holes which correspond to a dual superconducting phase. I concentrate on the case 
of black holes with a hyperbolic horizon because their properties can be understood analytically. 

In section|2]l discuss scalar, gravitational and electromagnetic perturbations of an AdS Schwarzschild 
black hole analytically calculating the QNM spectrum in the high frequency regime. In section [3] I 
calculate the QNM spectrum analytically in the low frequency regime and discuss its relevance to 
the hyrdodynamic behavior of the dual gauge theory fluid on the boundary. In sectional introduce 
hairy black holes and discuss their phase transition in the case of a hyperbolic horizon which can be 
understood analytically. The dual gauge theory corresponds to a superconductor whose properties 
can be calculated via electromagnetic perturbations. Finally, I conclude with section |5] 

2 Perturbations 

In this section I discuss scalar, gravitational and electromagnetic perturbations of an AdS Schwarzschild 
black hole in d dimensions analytically calculating the QNM spectrum in the high frequency regime. 
Low overtones will be discussed in the next section. 
The metric of an AdS Schwarzschild black hole is 

d^'' = -ii+K~^^dt^ + l^^!^ + ^d^h-2 (2) 



/2 ^/-3 / r^- ,1^ 2/i 

I shall choose units so that the AdS radius I = \. The horizon radius and Hawking temperature are, 
respectively. 



r 
The mass and entropy of the hole are, respectively. 



Anr4 



M = id-2){K + rl)^VolCLKj-2) , S = "^VolCLK^d-i) (4) 



The parameter K determines the curvature of the horizon and the boundary of AdS space. For K = 
0, +1 , — 1 we have, respectively, a flat (M"^^^), spherical (S'^^^) and hyperbolic (H^^^/F, topological 
black hole, where F is a discrete group of isometrics) horizon (boundary). 
The hai^monics on l^K,d-2 satisfy 

(V2 + A:2)T = (5) 

For A" = 0, /c is the momentum; for K = +1, the eigenvalues are quantized, 

k^ = l{l + d-?>)-5 (6) 

whereas for A' = — 1, 

where t, is discrete for non-trivial F. 5 = 0, 1 , 2 for scalar, vector, or tensor perturbations, respec- 
tively. 

According to the AdS/CFT correspondence, QNMs of AdS black holes are expected to corre- 
spond to perturbations of the dual Conformal Field Theory (CFT) on the boundary. The establish- 
ment of such a correspondence is hindered by difficulties in solving the wave equation governing 
the various types of perturbation. In three dimensions one obtains a hypergeometric equation which 
leads to explicit analytic expressions for the QNMs HI 13. In five dimensions one obtains a Heun 
equation and a derivation of analytic expressions for QNMs is no longer possible. On the other 
hand, numerical results exist in four, five and seven dimensions |[3J|4J|3. 

2.1 Scalar perturbations 

To find the asymptotic form of QNMs, we need to find an approximation to the wave equation 
valid in the high frequency regime. In three dimensions the resulting wave equation will be an 
exact equation (hypergeometric equation). In five dimensions, I shall turn the Heun equation into a 
hypergeometric equation which will lead to an analytic expression for the asymptotic form of QNM 
frequencies in agreement with numerical results. 

2.1.1 AdSs 

In three dimensions the wave equation for a massless scalar field is 

\dr {r' (l - ^) 5.<I>) - -±-,df^+^^d-^ = (8) 

Writing the wavefunction in the form 



r2 



^ = e'^o,t-px)^,^^^^ ^^^ (9) 



r2 
the wave function becomes 

y^{y-\)[ij-\yv'^ ^(a^y^^P^y{y-\yv = ^ (10) 

to be solved in the interval < j < 1, where 

a) = = , v = = . (11) 

2r+ AnTu 2r+ AkTh 
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For QNMs, we are interested in the solution 

ny)=y{'^-yy''2Fi{i + i{(b+p),\ + ii(b-py,2;y) (i2) 

which vanishes at the boundary (y — )■ 0). Near the horizon (y — )• 1), we obtain a mixture of ingoing 
and outgoing waves, 



»F-A+(l-3;)-'«+A_(l-3;) 



+ia, , _ n±2ia)) 



' " r{\±i{(b+p))r{\±i{d)-p)) 

Setting A_ = 0, we deduce the quasi-normal frequencies 

d) = ±p — in , « = 1,2,... (13) 

which form a discrete spectrum of complex frequencies with 3 <y < 0. 

2.1.2 AdSs 

Restricting attention to the case of a large black hole, the massless scalar wave equation reads 

l5.(rV(r) 5.0)- ^5,^0 -iv^<I> = 0, /(r) = 1 - ^ (14) 



Writing the solution in the form 



j2 
^ = e'^""-P-^'^^{y) , y = ^ (15) 
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the radial wave equation becomes 

(/-l)(K/-im'+(^/-^(/-l))^ = (16) 

For QNMs, we are interested in the analytic solution which vanishes at the boundary and behaves as 
an ingoing wave at the horizon. The wave equation contains an additional (unphysical) singularity 
at J = —1, at which the wavef unction behaves as *P '^ (y + 1)^®/^. Isolating the behavior of the 
wavefunction near the singularities y = ±1, 

'i'iy) = (y-ir^^'iy + lf^^'F^iy) (17) 

we shall obtain two sets of modes with the same 3(y, but opposite 9?(y. 
F± (y) satisfies the Heun equation 



v(/-l)f±" + {(3-^<i))j^-^<i>J-l}fi 



to be solved in a region in the complex j-plane containing jy] > 1 which includes the physical 
regime r> r+. 



For large cb, the constant terms in the polynomial coefficients of F' and F are small compared 
with the other terms, therefore they may be dropped. The wave equation may then be approximated 
by a hypergeometric equation 

(/-1W+{(3-^«').-^«'}44(±^=f1-.)f, = (19, 

in the asymptotic limit of large frequencies cb. The acceptable solution is 

Fo{x) = 2Fi{a+,a^;c;{y + \)/2) , a± = l-i^a)±l , c = |±^a> (20) 

For proper behavior at the boundary {y — ;■ oo), we demand that F be a polynomial, which leads to 
the condition 

a+ = —n , n=l,2,... (21) 

Indeed, it implies that F is a polynomial of order n, so as j — )• oo, F ~ 3^" ~ 3^^"+ and *F ~ 
y-ia/4y±m/4y-a^ ~ y-^, as expected. 

We deduce the quasi-normal frequencies 161 

(b = -^=2n{±l-i) (22) 

4nTH 

in agreement with numerical results. 

It is perhaps worth mentioning that these frequencies may also be deduced by a simple mon- 
odromy argument 161 . Considering the monodromies around the singularities, if the wavef unction 
has no singularities other than y = ±1, the contour around y = +1 may be unobstructedly deformed 
into the contour around y = —I, which yields 

^(l)^(-l) = 1 (23) 

Since the respective monodromies are ^(1) = e'^®/^ and ^(—1) = e^^'^^^^, using 3ft) < 0, we 
deduce cb = 2n{±l — i), in agreement with our result above. 

2.2 Gravitational perturbations 

Next I consider gravitational perturbations. For definiteness, I concentrate on the case of spherical 
black holes (K = +1). I shall derive analytic expressions for QNMs [7] including first-order cor- 
rections fH. The results are in good agreement with results of numerical analysis tOJ. Extension to 
other forms of the horizon is straightforward ITOl . 

The radial wave equation for gravitational perturbations in the black-hole background (O can 
be cast into a Schrodinger-like form, 

-^+V[rir,)]^ = CO^W, (24) 

in terms of the tortoise coordinate defined by 

'" ' (25) 



dr f{r) 

The potential V for the various types of perturbation has been found by Ishibashi and Kodama ITTI . 
For tensor, vector and scalar perturbations, one obtains, respectively. 



Vs(r) 



Vv{r)=f{r) 

fir 



+ d-3) , id-2){d-4)f{r) rf'ir) 



■ + 



4^2 



2{d-3) 



4r2 



e{e + d-3)-{d-2) + 



id-\)id-2)^ 



-2 



d{d-\Y{d-2yn^ 6{d-\){d-2y{d-4)[i{i + d-3)-{d-2)]n 



g2f.2d-i 



gljJ-S 



{d-A){d-€)[i{t + d-3)-{d-2)fr^ 2{d - if {d - 2^ ^^ 



/?2 



rid-9 



(27) 



+ 



A{d-\){d-2){2d'^-\\d+n)[l{l + d-3)-{d-2)\^'^ 

y2d—6 

{d-lf{d-2f{d-4){d-6)n^ 6{d-2){d-6)[£{i + d-3)-{d-2)fiJ. 



f.2d—6 



yd-?, 



6{d-l){d-2f{d-4)[i{i + d-3)-id-2)]^ 

fd—3 



+4[i{i + d-3)-{d-2)]^ + d{d- 2) [i{i + d-3)-{d-2){' 



Near the black hole singularity (r ~ 0), 

1 £/t 

Vj = -—^ H -n 

^n {-2[d-2)]x\—^ 



d-\ 

d-1 



{d-3f 



l + d-3) 



Vv 



M/ 



d-2 



'•''. [-2(</-2)fl]n 



+ ■•■■ ■■*'-^ = #7^ + ''^l?^2-^- <'» 

^ d^-Sd+l3 lU + d-3) 
Vs = -— + . ,,, . ,1 r, '^-^ + ..., (30) 

4(2J-5)(^-2) '^ {d-l){d-2Y ' ^^^^ 

I have included only the terms which contribute to the order I am interested in. The behavior of the 
potential near the origin may be summarized by 

/-I 



and 



where 



(2J3 - 24J2 + 94J _ 1 16) ( J2 _ 7 J + 14) [£(^ ^d-3)-{d-2)] 



V 



4ri 



+ £/n ''"' + 



d-l 
d-1 



(32) 



where 7 = (2) for scalar and tensor (vector) perturbations. 
On the other hand, near the boundary (large r), 



W 



ji-1 



4(r*-rJ 



r* 



dr 

'0 7w 



(33) 



where 7t>o = J — 1 , (i — 3 and J — 5 for tensor, vector and scalar perturbations, respectively. 
After rescahng the tortoise coordinate (z = G)r*), the wave equation to first order becomes 



J^ + w"?^J^ »F = 



(34) 



where 

^0 = ^ 



/-I 
■' 1 



jri = -i/z"3^. (35) 



4^2 

By treating ^ as a perturbation, one may expand the wave function 

»P(z)=»Po(z) + «-^»Pi(z) + ... (36) 

and solve the wave equation perturbatively. 
The zeroth-order wave equation, 

^o*I'o(z)=0, (37) 

may be solved in terms of Bessel functions, 

*I'o(z)=AiVi-/;(z)+A2V2^.(z). (38) 

2 2 

For large z, it behaves as 

^{){z) ~ A/-[Aicos(z-a+)+A2sin(z-a+)] 

- — "— (Ai - /A2)e-'"+e'^ + -^ (Ai + /A2)e+'«+e-'^ 



271 V27r 

where a± = f (1± j). 

At the boundary (r — ;■ oo), the wavefunction ought to vanish, therefore the acceptable solution is 



»Fo(r,) = B^(o{r^ - r*) J i^ («(''* " ^*)) (39) 

2 

Indeed, *F — )■ as r* — )■ f* , as desired. 
Asymptotically (large z), it behaves as 

»F(r,)~A/-Bcos[w(r*-r*)+j3] , j3 = ^(l+j.c) (40) 

V TT 4 

This ought to be matched to the asymptotic form of the wavefunction in the vicinity of the black- 
hole singularity along the Stokes hne 3z = 3 (wr* ) = 0. This leads to a constraint on the coefficients 
Ai, A2, 

Aitan(a)r*-j3-a+)-A2 = 0. (41) 

By imposing the boundary condition at the horizon 

»F(z)~e'^ , z^-00, (42) 

one obtains a second constraint. To find it, one needs to analytically continue the wavefunction 
near the black hole singularity (z = 0) to negative values of z. A rotation of z by —% corresponds 
to a rotation by — -^^ near the origin in the complex r-plane. Using the known behavior of Bessel 
functions 

/^(e-'^^z) = e-''^Vv(z) , f^vie-'^'z) = e'''''Ny{z) - 2/ cos nvJy{z) (43) 

for z < the wavefunction changes to 

»I/o(z) = e-'''^J+''^/^V^z { [Ai - i{\+e"'J)A2] J^{-z) +A2e''^^'A^. (-z)} (44) 



whose asymptotic behavior is given by 



-r7t(;+l)/2 -iTz(j+l)l2 

'i' ~ ^^ [Ai - i{\+le^''')A2] e-'"- + - — j=^ [Ay - ^^2] e'^ 



V271 "■ ' V2n 

Therefore one obtains a second constraint 

Ai-i{l+2eJ'")A2 = . 

The two constraints are compatible provided 

1 -i{\+2eJ'") 

tan(a)r* -j8-a+) -1 

which yields the quasi-normal frequencies Q 



cor, = -{2 + j + j^) -tan ^-^—^ + nK 



(45) 



(46) 



(47) 



(48) 



The first-order correction to the above asymptotic expression may be found by standard pertur- 
bation theory |^. To first order, the wave equation becomes 



The solution is 



'i'i{z) = VzNj{z) / dz' 



j^»Fi + jri»Fo = o 



V7Jj{z')Jn'i'o{z') rz ^/?Nj{z'Wl'i>o{z') 

' ^ ^zJiiz) dz' ^ 

2 JO 



w 



(49) 



(50) 



where W = 2/ 71 is the Wronskian. 

The wavefunction to first order reads 



^<{z) = {A,[\-b{z)] -A2a2{z)}V~zJi{z) + {A2[l+b{z)] + A,a,{z)} ^zN dz) (51) 



where 



«i(z) 

aiiz) 

biz) 



^ a,-^- dz! z! '-Uj{z')Jj{z) 
2 Jo 22 

-^(O-— / dz! z! '-^N.i{z!)Nj{z') 

2 Jo 2 2 



d-3 



^ (o ^- I'dz'z' "-Ujiz'jNjiz') 
2 Jo 22 



and £/ depends on the type of perturbation. 
Asymptotically, it behaves as 



*F(z) ~ Y - [A\ cos(z - a+) +A2sin(z - «+)] , 



where 



A[ = [I - b]Ai - a2A2 , A2 = [l+^]A2 + aiAi 
and I introduced the notation 



(52) 
(53) 



ai =01(00) , (22 = 02(00) 



(54) 



The first constraint is modified to 



Explicitly, 



Aitan(cof*-j8-a+)-A2 = (55) 



[(1 -b)tan{(on -j3-a+)-ai]Ai - [1 +Z7 + a2tan(G)r* -/3 -a+)]A2 = (56) 



To find the second constraint to first order, one needs to approach the horizon. This entails a rotation 
by —71 in the z-plane. Using 



a.ie-'^'z) 



a2{e-'\) 



z-'^'T^e-'''iax{z) 






e'''ia2{z) - 4cos2 ^a^ {z) - 2/(1 + e'''^)b{z) 



bie-'^'z) = e-'''—2[biz)-i{\+e-'''J)aiiz)] 
in the limit z — > — oo one obtains 

^{z) ~ -ie-'J^'/^Bi cos(-z - a+) - /e'^'^^/^g^ sin(-z - a+) 



(57) 



where 



Bi = Ai-Aie-'''^[b-i{l+e-'''J)ai] 



-A2e" 



e+'^'Jai - 4cos2 ^aj - 2/(1 + e+''^-'> 



B9 






Therefore the second constraint to first order reads 

[1 - e-'^^ef (/^j + ^Ai - [/(I + le'^'J) + e-'^^^ef ((1 + e'''J)di + e''^-'a2 - /^)]A2 = 
Compatibility of the two first-order constraints yields 

1+^ + ^2 tan(wr* - j3 - a+) /(I + le^^J) + e-'""^ ((1 + e''^-'>i + e''^-'a2 - ib) 
(1 -l7)tan(G)r* -/3-a+)-ai 1 -e^'^d^{idi +b) 

leading to the first-order expression for quasi-normal frequencies, 

^ /^ . . \ 1 , ^ 

(Or^ = —{2 + j + joo) + —ln2 + nn 



(58) 



(59) 



where 



-i{«' 



;„ d—3 — ;„ d—3 ;„ d—3 

Gib -lie '^^r^b-9d[+e "^~ai+a2-e 3^a2 



ai 

d2 
b 



nil 

4 V2^ 



rf-3 
rf-2 



^(^-2)^(2 + 2(d-2)> 

r2( d-i \Y(J I ^-1 - 

^ V2(rf-2)^^ V2 ^ 2(d-2)' 



, ^ 7r(<i-3) K f . d-3 
l+2cot:r^^ — TT-cot- I -j + 



2{d-2) 2 



d-2 



a\ 



nid-3) 



rf-3 , 



Thus the first-order correction is ~ ^(n^^?^). 

The above analytic results are in good agreement with numerical results ^ (see ref. HI for a 
detailed comparison). 

2.3 Electromagnetic perturbations 

The electromagnetic potential in four dimensions is 

Vem = ^^^/W. (60) 

Near the origin, 

where 7 = 1. Therefore a vanishing potential to zeroth order is obtained. To calculate the QNM 
spectrum one needs to include first-order corrections from the outset. Working as with gravitational 
perturbations, one obtains the QNMs 

(Or^ = nn -- Inn + — In {2{l + i)£/^/f^) , js^ = ^^—==L (62) 

Notice that the first-order correction behaves as Inn, a fact which may be associated with gauge 
invariance. 

As with gravitational perturbations, the above analytic results are in good agreement with nu- 
merical results 19J (see ref. L8J for a detailed comparison). 



3 Hydrodynamics 



There is a correspondence between .yV = 4 Super Yang -Mills (SYM) theory in the large N limit 
and type-IIB string theory in AdSs x S^ (AdS/CFT correspondence). In the low energy limit, string 
theory is reduced to classical supergravity and the AdS/CFT correspondence allows one to calcu- 
late all gauge field-theory correlation functions in the strong coupling limit leading to non-trivial 
predictions on the behavior of gauge theory fluids. For example, the entropy of o/K = 4 SYM theory 
in the limit of large 't Hooft coupling is precisely 3/4 its value in the zero coupling limit. 

The long-distance (low-frequency) behavior of any interacting theory at finite temperature must 
be described by fluid mechanics (hydrodynamics). This leads to a universality in physical proper- 
ties because hydrodynamics implies very precise constraints on correlation functions of conserved 
currents and the stress-energy tensor. Their correlators are fixed once a few transport coefficients 
are known. 

3.1 Vector perturbations 

I start with vector perturbations and work in the li-dimensional Schwarzschild background ^ with 
K = +\ (spherical horizon and boundary). It is convenient to introduce the coordinate [il2il 



(^) 



10 



(63) 



The wave equation becomes 



rf-4 /^. . / rf-4 /^ 



(0 



{d - 3yu—^f{u) u—^f{u)Y + Vs/{u)^ = a)^^> , (£) = — 



where prime denotes differentiation with respect to u and I have defined 



-/ N fir) 2 / 1 — M 

/(m) = ^ = 1-m— (m 



r± 



where L^ = -^^-^ — - . 



., {d-2){d-4) 2 . , {d-\){d-2)h+Jr 



(64) 



(65) 



(66) 



First I consider the large black hole limit r_|_ — )■ oo keeping <y and L fixed (small). Factoring out 
the behavior at the horizon (m = 1) 



the wave equation simplifies to 



where 



^'=(\-u)-'^F(u) 



s^F" + .^^F' + '^^lF=0 



2rf-8 , d-l . 



-{d-3yu^^{\-u^ 



d-l , d-5 



'^(bl 



-{d-3)[d-A-{2d-5)u—]u—-2{d-3) 

., , {d-2)[d-A-3{d-2)ui^^] _2 






Id-i d-l , 

. ico W-i (1 -Ud-i ^ 

d—\ \—u 



1 -Ud-i 



— + {d-3y 



{d-l) 



l-u) 



i(0 [d - 4 - {2d - 5)ud-i\u-'-i ^ ^^ U''-i{\-ud-\ 



d — I \ — u 

One may solve this equation perturbatively by separating 

{M)+-m)F=o 

where 

^qF = £/F" + ^QF' + 'ifo,oF 

J^F = (^d,-^())^' + C^d,X-^o,o)F 

Expanding the wavefunction perturbatively, 

F = Fo + Fi + ... 



d-l {l-uf 



(67) 



(68) 



(69) 



(70) 
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at zeroth order the wave equation reads 

^0^0 = (71) 

whose acceptable solution is 



Fo = m2(''-3) (72) 

being regular at both the horizon (u = 1) and the boundary (m = 0, or *F '^ r~ t~ — ;• as r — )■ oo). 
The Wronskian is 

^= .-4, ^ .-, (73) 



and another linearly independent solution is 



/T 



Fo = Fo -J (74) 

which is unacceptable because it diverges at both the horizon (Fq ^ ln(l — u) for m ss 1) and the 

rf-4 rf-4 

boundary (Fq ~ u ^'■''-^'i for m « 0, or *F ~ r 2 — )■ 00 as r — )• 00). 
At first order the wave equation reads 

JToFi = -J^Fo (75) 

whose solution may be written as 

F? J £/W 



Fi=Fol ^i ^^Z7^ (76) 



The limits of the inner integral may be adjusted at will because this amounts to adding an arbitrary 
amount of the unacceptable solution. To ensure regularity at the horizon, choose one of the limits 
of integration at m = 1 rendering the integrand regular at the horizon. Then at the boundary {u = 0), 

. /■iFoJ^fo 
F[=Fo ...^ + regular terms (77) 

The coefficient of the singularity ought to vanish, 

■'^=« ™ 

siW 
which yields a constraint on the parameters (dispersion relation) 

aoL^-miO)- 3201^ = (79) 



After some algebra, one arrives at 



J-3 
ao = - — - , ai = J - 3 (80) 

d—\ 



The coefficient a2 may also be found explicitly for each dimension d, but it cannot be written as a 
function of d in closed form. It does not contribute to the dispersion relation at lowest order. E.g., 
for J = 4,5, one obtains, respectively 
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Eq. ( |79l ) is quadratic in a) and has two solutions, 



COO 



Wl 



^-3 n 



(i — 1 " 82 

In terms of the frequency (O and the quantum number (., 



m + d-3) (Oi d-3 m + d-3) 

{d-l)r+ r+ 32 {d-l)ri 



(82) 



(83) 



The smaller of the two, coq, is inversely proportional to the radius of the horizon and is not included 
in the asymptotic spectrum. The other solution, COi , is a crude estimate of the first overtone in the 
asymptotic spectrum, nevertheless it shares two important features with the asymptotic spectrum: 
it is proportional to r+ and its dependence on i is ^(1/r^). The approximation may be improved 
by including higher-order terms. This increases the degree of the polynomial in the dispersion 
relation ( |79l ) whose roots then yield approximate values of more QNMs. This method reproduces 
the asymptotic spectrum derived earlier albeit not in an efficient way. 

To include finite size effects, I shall use perturbation theory (assuming l/r+ is small) and replace 
=^iby 

jif' = j^i + -^jr+ (84) 

where 

J^+F = £^+F" + ^+F' + '^+F (85) 

The coefficients may be easily deduced by collecting ^(1/r^) terms in the exact wave equation. 
One obtains 



-2{d-3fu^{l-u) 
-{d — 3)u 



{d - 3)(2 - 3m) -{d- 1) ^M^J 



1 -W-s 



^+ 



d-2 



\ — U d-\ 

d-4-(2d-5)u-{d- 1) t^u^ 



1 -W-a 



Interestingly, the zeroth order wavefunction Fq is an eigenfunction of J^, 

J^+Fo = -{d-2)Fo 
therefore the first-order finite-size effect is a simple shift of the angular momentum operator 

d-2 



L2 



P- 



(86) 



(87) 



The QNMs of lowest frequency are modified to 

.e{£ + d-3)-{d-2) 



(Oq 



For fif = 4,5, we have respectively. 



(Oo 



id-\)r+ 



3r^ 



+ffil/rl) 



4r4 



(88) 



(89) 
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in agreement with numerical results |[9l[T3l. 

One deduces from (1881 ) the maximum lifetime of the vector modes, 



471 

Vax = -r7)/ (90) 

a 



In the case of a flat horizon {K = 0), 



(OQ = -i ,, \^ (91) 

[d-\)r+ 

which leads to the diffusion constant 

In the case of a hyperbolic horizon (K = — 1), a similar calculation yields ifTOll 

((3 — l)r+ |G)o| ((3 — 1)^ 

It follows that for J = 5, these modes live longer than their spherical counterparts which is important 
for plasma behavior. 

3.2 Scalar perturbations 

Next I consider scalar perturbations which are calculationally more involved but phenomenolog- 
ically more important because their spectrum contains the lowest frequencies and therefore the 
longest living modes. For a scalar perturbation we ought to replace the potential Vy by 



t>s(") 



m+ { 1 + — I M 

3 



+{d-4){d-6)m^u-— + {d-2)^il + —] u^ 

+2{2d^-nd + n)m(\ + \\ u^ 

{d-A){d-€)(\ + ^)^ f 1 ^ 

+ 2"^ ^i^M^ - 3((/ - 2){d - G)n? 1 + — ) M 

6(J-2)(J-4)m(l + 4) ^2 
-2 ^ ^u + 2{d - \){d -2)m^ +d{d -2) — 



(94) 



where m - n ^(M-3)-{d-2) _ 2{e+d-2){i-l) 
Where m - 2 (^„i)(^_2)ri " (J-i)(rf-2)4 " 
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In the large black hole limit r+ — ;• oo with tn fixed (small), the potential simplifies to 



f/(0). ^ l-U^ 



VPiu) 



d{d-2)u'^^ -6{d-2){d-A)mu^^ 



+{d-A){d- 6)m^u-^ + id- ifw" 

+2{2d^ -nd + \S)mu'^ - 3{d -2){d- 6)m^u + 2{d -\){d- 2)nr' 



(95) 



The wave equation has an additional singularity due to the double pole of the scalar potential at 
u = —tn. It is desirable to factor out the behavior not only at the horizon, but also at the boundary 
and the pole of the scalar potential. 



Then the wave equation reads 
where 



»P =!-„-'— _ F(u) 

^ ' m + u ^ ' 



.s^F" + ,^s,F' + %j,F = 



(96) 
(97) 



-{d-3yu^r^{\-u^) 
-[d-djW-i (l-M^-s) 



d-4 2{d-3) 



m + u 



2d-» d-1 , 

, d-l, d-5 , ,9 im M ''-S (1 — M<'-3 

-(d-3){d-A-{2d-5)u—^\u—^-2{d-3f ^ 



%] 



2ri-8 rf-1 

-M''-3 (l-W-S 



d—\ \—u 

{d-2){d-4) {d-3){d-4) 2{d-3)'- 



+ ■ 



4u^ u{m + u) (m + uY 

2d-S ±2ls ' 

d-4-{2d-5)u— \u—+2{d-3)- — 



(i-4 d-3 



2u m + u 



JO) \d - 4 - (2d - 5)u'i-^\u'i-^ ^ '^ m ^-s (1 - m^-s 

-(^(3 — jj— - - - [a — 5) 



+ - 



d-\ 

vf\u\-(iY 



\ — u 



d-\ {\-uf 



d-\ 
1 -Ui-^ 



+ (d-3f 



or w-i (1 — w 



3 



(J- 1)2 (l-i<)2 

I shall define the zeroth-order wave equation as J^)Fo = 0, where 

J^oF = ^F" + ^qF' 

The acceptable zeroth-order solution is 

FnfM) = 1 



(98) 



(99) 



which is plainly regular at all singular points (m = 0, 1,— m). It corresponds to a wavefunction 

d-A 

vanishing at the boundary (^ r^ r 2 as r — )■ 00). 
The Wronskian is 



W-- 



im + u) 



2d-8 . d-\ , 

W-S (1 -M^-S 



(100) 
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and an naccep table solution is Fq = JW. It can be written in terms of hypergeometric functions. 

^- rf-5 rf-6 

For d > 6, it has a singularity at the boundary, Fq ^ u ^-3 for m ss 0, or ^F ~ r 2 — >. 00 as r — )• 00. 
For J = 5, the acceptable wavef unction behaves as r^'/^ whereas the unacceptable one behaves as 
r^^'^lnr. For d = A, the roles of Fq and Fq are reversed, however the results still valid because the 
correct boundary condition at the boundary is a Robin boundary condition |[T2llT4]| . Finally, note 
that Fq is also singular (logarithmically) at the horizon {u = 1). 

Working as in the case of vector modes, one arrives at the first-order constraint 



' % 



'" (101) 



/O £/W 
because J^\Fq = {■3Sq, — ^o)Fq + 'j^aFo = ^^. This leads to the dispersion relation 



ao-ai/a)-a2ai^ = (102) 



After some algebra, one obtains 



d-\ l + (d-2)m d-3 1 , 

«o = -5 nV M2 ' ai= , a2 = -{l + (9m} (103) 

For small th, the quadratic equation has solutions 



% - -^—T- "J ± A/ —T- ^ (104) 



^+ .J — 3^, d—\^ 



related to each other by G)q = — G)q , which is a general symmetry of the spectrum. 

Finite size effects at first order amount to a shift of the coefficient ao in the dispersion relation 

ao-^ao + ^a+ (105) 

After some tedious but straightforward algebra, we obtain 

a+ = 4{l + 0(m)} (106) 

m 



The modified dispersion relation yields the modes 



.±., .^-3 _ ld-\ ^ 



In terms of the quantum number £, 



,^ ,^+d-^y^d-2;i l i{i+d- 3) 

'« ~"'^^"^^ (J-l)(J-2)r+ ^y^d^ 



co^K-i{d-3) ^-j — / ; ^ifcW ^ ' ^ (108) 



in agreement with numerical results |[T3l . 

Notice that the imaginary part is inversely proportional to r+, as in vector case. In the scalar 
case, we also obtained a finite real part independent of r+. 

The maximum lifetime of a gravitational scalar mode is found from (11081 ) to be 

d-2 

AtzTh (109) 



{d-3)d 
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In the case of a flat horizon (K = 0), one obtains 



k d-3 



2 



(o = ±—^=-i— — -— — - — r (110) 



"V^32 {d-\){d-l)r+ 
showing that the speed of sound is 



' = 7ki <"" 



as expected for a CFT and the diffusion constant is 



d-3 1 



For a hyperbolic horizon {K = — 1), a similar calculation yields ifTOl 



In the physically relevant case d = 5, evidently the K = —I scalar modes live longer than any other 
modes, which is important for plasma behavior. 

3.3 Tensor perturbations 

Finally, for completeness I discuss the case of tensor perturbations. Unlike the other two cases of 
gravitational perturbations, the asymptotic spectrum of tensor perturbations is the entire spectrum. 
To see this, note that in the large black hole limit, the wave equation reads 

-{d- 3f{u^ - u^)'¥" -{d-3) [{d - 4)u^ - {2d - 5)u^]'¥' 
(f2 d{d-2) _^ {d-lf (iP- \„, 

I 4 4 l-u—i) 

For the zeroth-order equation, we may set L = = O). The resulting equation may be solved exactly. 
Two hnearly independent solutions are {^ = Fq at zeroth order) 

Fo{u)=u^(d^) , Fo(m) = M~2(3^1n('l-M3^') (114) 

Neither behaves nicely at both ends (u = 0, 1). Therefore both are unacceptable which makes it 
impossible to build a perturbation theory to calculate small frequencies which are inversely propor- 
tional to ro. This negative result is in agreement with numerical results |l9l[T3l and in accordance 
with the AdS/CFT correspondence. Indeed, there is no ansatz that can be built from tensor spherical 
harmonics T;y satisfying the linearized hydrodynamic equations, because of the conservation and 
tracelessness properties of Tjj. 

3.4 Hydrodynamics on the AdS boundary 

The above results in the bulk dictate the hydrodynamic behavior of the dual gauge theory fluid on 
the conformal boundary. To see the correspondence, one needs to understand the hydrodynamics in 
the linearized regime of a (i — 1 dimensional fluid with dissipative effects. The fluid lives on a space 
with metric 

dsl = -dt^ + dLlj_2 (115) 
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The hydrodynamic equations are simply the requirement that the stress-energy momentum tensor 
be conserved, 

V^T^'' = (116) 

As the duality corresponds to a conformal field theory one must also demand scale invariance which 
implies 

r^=0, e = {d-2)p , C=0 (117) 

where e, p and ^ are the energy density, pressure and bulk viscosity of the fluid. In the rest frame of 
the fluid, the velocity field is m^ = ( 1 , 0, 0, 0) and the pressure po is constant. Consider a perturbation 



[l,u') , p = pQ + 5p 



Applying the hydrodynamic equations, one obtains 



{d- l)pod,u' + d'5p-ri 



{d-2)d,5p + {d-l)poViu' 
VJVju'+K{d - 3)m' + ^^d'iVjuJ) 








where I used the curvature tensor Rij = K{d — 3)gij. 
For vector perturbations, consider the ansatz 

dp = 0, u' = ^ve''""T 

where V is a vector harmonic. 

The hydrodynamic equations imply 

-i(o{d-l)po + ri [4-^(^-3)] =0 

Using 

R. = (d-2)^— = ^^^ ''^ 

po sM s K + rl 

with (O from the gravity dual, one obtains for large r+, 

1 = 1- 

s Ak 

which is the standard value of the ratio in gauge theory fluids with a gravity dual ifTSll . 
For scalar perturbations, consider the ansatz 

u' = ^se-""d'S , dp = mse-'""§ 

where S is a scalar harmonic. 

The hydrodynamic equations imply the system of equations 

(d - 2)i(0^s + {d- l)pokj.<2/s = 



The determinant must vanish, 

{d - 2)i(0 



-i(o{d -l)po- 2{d - 3)Kri + 2riks 



{d-\)pokj 



2^-3 



d-2 







1 -i(o{d-l)po-2{d-3)Kri+2rikj^ 

Arguing along the same lines as for vector perturbations, we arrive at 

1 = 1- 

s An 

which is the same result as the one obtained with vector QNMs. 



(118) 



(119) 

(120) 

(121) 
(122) 

(123) 
(124) 



(125) 



(126) 



(127) 
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3.5 Conformal soliton flow 

The above results have been applied to the study of the quark-gluon plasma which forms in heavy 
ion collisions (at the Relativistic Heavy Ion Collider (RHIC) and elsewhere). In the case of a 
spherical horizon (K = -\-l), the boundary of spacetime is S^ x M. This may be conformally mapped 
onto a flat Minkowski space. Then by holographic renormalization, the AdSs-Schwarzschild black 
hole is dual to a spherical shell of plasma on the four-dimensional Minkowski space which first 
contracts and then expands (conformal soliton flow) ifTSll . 

Quasi-normal modes govern the properties of this plasma with long-lived modes (i.e., of small 
3 ft)) having the most influence. For example, one obtains the ratio 

V2 1 (0^-40(0^ + 12 . K(0 
671 0)3-40) 2 

where V2 = (cos 20) evaluated at 6 = ^ (mid-rapidity) and averaged with respect to the energy 
density at late times; d = ['[f/'.il is the eccentricity at time t = 0. Numerically, ^ = 0.37, which 
compares well with the result from RHIC data, ^ ^ 0.323 1 16|. 

Another observable is the thermalization time which is found to be 

^ = Tk!^ ^ ^^^ ^ 0.08 fm/c , rpeak = 300MeV (129) 

Z|JiO)| o.oypeak 

not in agreement with the RHIC result T ~ 0.6 fm/c |fT7l|, but still encouragingly small. For com- 
parison, the corresponding result from perturbative QCD is T > 2.5 fm/c lfT8l[T9ll . 

In the case of a hyperbolic horizon (topological black hole; ^ = — 1), one needs to work with a 
conformal map from H^^^/F x M to a (J — 1)— dimensional Minkowski space. Finding an explicit 
form of this map for J = 5 involves a considerable amount of numerical work. However, it is 
important that one consider this case because the modes of hyperbolic black holes live the longest 

Eol. 



4 Phase transitions 

In this section I discuss hairy black holes in asymptotically AdS space and their duals. At low 
temperatures, an instability leads to symmetry breaking and the formation of a dual superconductor. 
ELectromagnetic perturbations of the black hole determine the conductivity in the bulk. First I 
review the case of a flat horizon {K = 0) f20| and then I discuss the case of hyperbolic horizon 
{K = —I) where exact analytical results are obtained 1.211 . 

4.1 K = 

Consider a scalar *P of mass m^ = —2, which is above the Breitenlohner-Freedman (BF) bound 
coupled to an electromagnetic potential A^ in 3-1-1 dimensions. The Lagrangian density is 

^ = _l5^»i/5^»j/ + »i/2 _ ^F^^F^" - y^^(a^0 -A^){d^d-A^') (130) 

where is a Stiickelberg field, q is an arbitrary parameter which can be throught of as the electric 
charge of the scalar field *P (one may instead turn ^ into a complex scalar field of charge q coupled 
to an electromagnetic potential in a standard fashion). 
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The Lagrangian density ( 11301 ) is invariant under the U{1) gauge transformation 

A^^A^ + d^co , e^e + (0 (131) 

To fix the gauge, set 

0=0 (132) 

Working in the probe limit {q — ;■ oo) in which there is no back reaction to the metric, assume that the 
fields propagate in the black hole background (O with d = A and K = 0. The radius of the horizon 
and Hawking temperature are, respectively, 

r+ = (2M)^/^ r = ^ (133) 

Assuming spherical symmetry and an electrostatic potential Aq = <I>(r), the field equations yield 
two coupled non-linear differential equations ll20l 



*!'"+( 4 + - 1*1''+ (^ I *I'+^*I' = 



/ rj \fj f 



where I set ^ = 1 and 



2 2*F^ 
O" + -<!)' O = (134) 

r f 



f{r)=r'-^ (135) 

r 



As r — )• oo, one obtains the boundary behavior 

»P = + ^ + ... , <I> = <I>(0) + + ... (136) 

r r'^ r 

where one of the *F(') =0 (/ = 1,2) for stability, O^"^) is the chemical potential and O^i) = — p 
(charge density). 

Below a critical temperature Tq a condensate forms, 

(^,.) = V2»F(') (137) 

of an operator of dimension A = /. 

At r = To, one may set *F = in the equation for <I> and deduce (using <I>(r+) = 0) 

^ = p(^--\ (138) 

Then the equation for ^ turns into an eigenvalue problem which yields 

To p:i 0.226 Vp , 0.118VP 

depending on the boundary conditions. 

To study the properties of the dual CFT, apply an electromagnetic perturbation. It obeys the 
wave equation 
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to be solved subject to the boundary conditions that it be ingoing at the horizon, A r^ f '(^/i^^T)^ 
and at the boundary (r — ;■ oo), 

4(1) 

A = a(°) + + ... (140) 

r 

Ohm's law yields the conductivity 

For T > Tq, *P = 0, therefore A ~ e""''" where r^ = Jdr/f{r) is the tortoise coordinate. It follows 
that 

a(w) = 1 (142) 

At low r, for (^i) / 0, we have 

Since r+ -^ 0, we obtain A ~ e'"''''*, where co' = \J (iP- - (^i)^- Therefore, for « < (i^i), 9ta = 0, 
i.e., we obtain a superconductor with a gap. 

4.2 ^ = -1 

Turning to the case of a hyperbolic horizon 121], choose a scalar *P of mass m^ = —2, as before, but 
conformally coupled with potential 

V(^) = ^^^ 
The system has an exact solution (MTZ black hole [22]) 

,2 / \ ^ 

ds^ = -fMTz{r)dt'^-^-—^?dG' , fMTZ = r'-[\-r-\ , (143) 

/Mrz('') V '" / 

with 



AnG r + rQ 



The temperature, entropy and mass are, respectively. 



1 / 1\ „ o ,^ ,, ., or+ 



T = -[r,--y W = ^(2r,-1), Mm.z = ^(..-1). (145) 

and the law of thermodynamics dM = TdS holds. 

At M = 0, the MTZ black hole coincides with the topological black hole with no hair (eq. Q 
with d = 4, K = -\), 



df- 
r2-l 
and an enhanced scaling symmetry (pure AdS space) emerges at the critical temperature 



^*Ls = -{r^ - l)dt'^ + -2 — 7 + r^dl? ( 146) 



7b = i^ (147) 

At this point there is a phase transition which can be seen by calculating the difference in free 
energies, 

^F=FTBH-FMTZ = -^7lH\T-TQf + ... , (148) 
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q/VG 


1 


3 


5 


5 


1.025 ±0.007 


1.52 ±0.03 


1.78 ±0.03 



Table 1 : The exponent d characterizing the low-temperature dependence of the superfluid density ris 

showing that there is a third-order phase transition at Tq. 

Perturbative stability of the MTZ black hole has also been demonsti^ated for T < To (M < 0) 
ll2n . Comparing with the flat case, note that here both ^F'^' and ^F^^' are non- vanishing, yet the 
MTZ black hole is stable. However this is true only if the mass is negative which is never the case 
with a flat horizon. Also, here the condensation of the scalar field has a geometrical origin and is 
due entirely to its coupling to gravity. 

Moreover, the heat capacities in the normal and superconducting (corresponding to the MTZ 
black hole) phases, respectively, as 7 — )• exhibit a power-law behavior 



Cn 



KG 

3V3G 



C^ — T 



(149) 



Since both ^^^' and *F(^' are non-vanishing, we have a multi-trace deformation of the CFT 123 
with a condensate 



37r3 
2G 



%^-T^) 



(150) 



It should be noted that the deformation does not break the global U{1) symmetiy because *P is a 
real field (see eq. (11301 )). 

To study the conductivity, apply an electromagnetic perturbation. It obeys the wave equation 
(II39I) which may be solved using first-order perturbation theory in q^, 



A = e-""-' + 



jmn 



lico 



/ui2fj\„-2i(ar. 



drW{r)e 



2im 



i\Jl2(J\ 



dr^>\r 



r+ 



(151) 



The conductivity to first order in q^ is 



a {CO) = ^-- = 1 - — / ^r»F2(r)^-2tor. 



(152) 



The superfluid density is found from 



91[<t(w)] ~ 7r«,.5(G)) , 3[a(G))]~ — , co^O. 



One obtains 



(153) 



ns = q^ I dr^>\r)- ^^ 



'■+ 



AnG r+ + ro 



a{TQ-Ty , a 



3nq^ 
~4G' 



(154) 



Near T = 0, 



«,(0)-«,(^)«-^^ 5 = 1 



n 



(155) 



In table [H this analytic prediction is compared against exact numerical results for various values of 
the charge q. Naturally, the agreement is best at small values of q. 
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q/VG 


?numerical 


Tanalytical 


O^numerical 


O^analytical 


0.1 


0.0020 


0.0024 


0.0225 


0.024 


0.5 


0.0538 


0.0597 


0.552 


0.589 


1.0 


0.187 


0.239 


2.196 


2.356 


2.0 


0.684 


0.955 


8.678 


9.425 


3.0 


1.325 


2.15 


20.35 


21.21 


5.0 


2.522 


5.97 


52.90 


58.90 



Table 2: Numerical vs analytical results for the normal and superfluid densities. 





Figure 1: The real part of the conductivity vs G) for qj \fG = 2 (left) and q/Vc 
0.0032,0.032,0.064. The lowest curve corresponds to the lowest temperature. 



5 (right) and T 



The normal, non-superconducting, component of the DC conductivity is 



Therefore, 



At low T, 



In 



fin = lim 9?[a(w)] , 



r'OO 

nn=lq^ / dr^>^{r)r, 

Jrj. 



"•+ 



tin^T^, y- 



AkG 



(156) 



(157) 



(158) 



This analytic result and the prediction for the parameter a determining the critical behavior of the 
superfluid density are compared against exact numerical results in table [2] Again, the agreement is 
best at small q. 

Figures [T] and [2] show the frequency dependence of the real and imaginary, respectively, parts 
of the conductivity. The real part of the conductivity becomes smaller as we increase the charge q. 
Unfortunately, numerical instabilities also increase and we have not been able to produce reliable 
numerical results above ql\fG = 5. The superconductor appears to be gapless. However, a gap is 
likely to develop above a certain value of the charge q, as indicated by the trend in the graphs as q 
increases. 
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Figure 2: The imaginary part of the conductivity multiplied hy (Ovs (O for q/y/G = 2 (left) and q/y/G -■ 
5 (right) and T = 0.0032,0.032,0.064. The uppermost curve corresponds to the lowest temperature. 

5 Conclusion 

The quasi-normal modes that govern perturbations of black holes in asymptotically AdS space 
are a powerful tool in understanding the hydrodynamic behavior of a gauge theory fluid at strong 
coupling. Here I focused on the analytic calculation of QNMs. I discussed both high overtones 
and low frequencies. I applied the results on gravitational perturbations to the understanding of the 
quark-gluon plasma produced in heavy ion collisions at RHIC and the LHC. I also considered hairy 
black holes whose electromagnetic perturbations allow one to analyze the conductivity of the dual 
conformal field theory and the phase transition to a superconducting state. I reviewed the case of 
a flat horizon and compared the results with those from black holes with hyperbolic horizon for 
which exact hairy solutions have been constructed (MTZ black holes |22]). In all these cases, only 
the low-lying QNMs were needed. It is unclear what physical role high overtones play. 



References 

[I] Cardoso, V., Lemos, J. P. S.: Scalar, electromagnetic and Weyl perturbations of BTZ black 
holes: quasi normal modes. Phys. Rev. D 63, 124015 (2001) 

[2] Birmingham, D., Sachs, I., Solodukhin, S. N.: Conformal Field Theory Interpretation of 
Black Hole Quasi-normal Modes. Phys. Rev. Lett. 88, 151301 (2002) 

[3] Horowitz, G. T, Hubeny, V. E.: Quasinormal Modes of AdS Black Holes and the Approach 
to Thermal Equilibrium. Phys. Rev. D 62, 024027 (2000) 

[4] Starinets, A. O.: Quasinormal Modes of Near Extremal Black Branes. Phys. Rev. D 66, 
124013 (2002) 

[5] Konoplya, R. A.: On quasinormal modes of small Schwarzschild-Anti-de-Sitter black hole. 
Phys. Rev. D 66, 044009 (2002) 

[6] Musiri, S., Siopsis, G.: Asymptotic form of quasi-normal modes of large AdS black holes. 
Phys. Lett. B 576, 309 (2003) 

[7] Natario, J., Schiappa, R.: On the Classification of Asymptotic Quasi-normal Frequencies for 
<i-Dimensional Black Holes and Quantum Gravity Adv. Theor. Math. Phys. 8, 1001 (2004) 

[8] Musiri, S., Ness, S., Siopsis, G.: Perturbative calculation of quasi-normal modes of AdS 
Schwarzschild black holes. Phys. Rev. D 73, 064001 (2006) 



24 



[9] Cardoso, V., Konoplya, R. A., Lemos, J. P. S.: Quasi-normal frequencies of Schwarzschild 
black holes in AdS space-times: A complete study on the asymptotic behavior. Phys. Rev. D 
68, 044024 (2003) 

[10] Alsup, J., Siopsis, G.: Low-lying quasinormal modes of topological AdS black holes and 
hydrodynamics. Phys. Rev. D 78, 086001 (2008) 

[11] Ishibashi, A., Kodama, H.: A Master Equation for Gravitational Perturbations of Maximally 
Symmetric Black Holes in Higher Dimensions. Prog. Theor. Phys. 110 701 (2003) 

[12] Siopsis, G.: Low frequency quasi-normal modes of AdS black holes. JHEP 0705, 042 (2007) 

[13] Friess, J. J., Gubser, S. S., Michalogiorgakis, G., Pufu, S. S.: Expanding plasmas and quasi- 
normal modes of anti-de Sitter black holes. JHEP 0704, 080 (2007) 

[14] Michalogiorgakis, G., Pufu, S. S.: Low-lying gravitational modes in the scalar sector of the 
global AdS4 black hole. JHEP 0702, 023 (2007) 

[15] Policastro, G., Son, D. T., Starinets, A. O.: From AdS/CFT correspondence to hydrodynam- 
ics. JHEP 0209, 043 (2002) 

[16] PHENIX Collaboration, Adare, A., et at: Scaling properties of azimuthal anisotropy in 
Au-i-Au and Cu-i-Cu collisions at ^/smn = 200 GeV. arXiv:nucl-ex/0608033 

[17] Arnold, P., Lenaghan, J., Moore, G. D., Yaffe, L. G.: Apparent Thermalization due to Plasma 
Instabilities in Quark-Gluon Plasma. Phys. Rev. Lett. 94, 072302 (2005) 

[18] Baier, R., Mueller, A. H., Schiff, D., Son, D. T.: "Bottom-up" thermalization in heavy ion 
coUisions. Phys. Lett. B 502, 51 (2001) 

[19] Molnar, D., Gyulassy, M.: Saturation of elliptic flow at RHIC: Results from the covariant 
elastic parton cascade model MPC. Nucl. Phys. A 697, 495 (2002) 

[20] Hartnoll, S. A., Herzog, C. P., Horowitz, G. T.: Building an AdS/CFT superconductor. Phys. 
Rev. Lett. 101, 031601 (2008) 

[21] Koutsoumbas, G., Papantonopoulos, E., Siopsis, G.: Exact Gravity Dual of a Gapless Super- 
conductor. JHEP 0907, 026 (2009) 

[22] Martinez, C, Troncoso, R., Zanelli, J.: Exact black hole solution with a minimally coupled 
scalar field. Phys. Rev. D 70, 084035 (2004) 

[23] Hertog, T., Maeda, K.: Black Holes with Scalar Hair and Asymptotics in N=8 Supergravity. 
JHEP 0407, 051 (2004) 



25 



